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Hamiltonians 



On 
O 
O 
(N 

> 

O 



43 ■ 
i ' 

cr 



> 
m 

(N 
O 



On 
O 



Boyan T. Torosov 1 and Nikolay V. Vitanov 1,2 

'Department of Physics, Sofia University, James Bourchier 5 blvd, 1164 Sofia, Bulgaria 
i Institute of Solid State Physics, Bulgarian Academy of Sciences, Tsarigradsko chaussee 72, 1784 Sofia, Bulgaria 

(Dated: October 6, 2009) 

We propose a technique for design of quantum Fourier transforms, and ensuing quantum algo- 
rithms, in a single interaction step by engineered Hamiltonians of circulant symmetry. The method 
uses adiabatic evolution and is robust against fluctuations of the interaction parameters as long as 
the Hamiltonian retains a circulant symmetry. 
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I. INTRODUCTION 

Quantum information processing is built upon se- 
quences of special unitary transformations. One of the 
most important of these is the quantum (discrete) Fourier 
transform (QFT), which is a key ingredient of many 
quantum algorithms [H, [3,13 , including Shor's factoriza- 
tion [J], the algorithms of Deutsch |B[ and Simon Q, 
order finding, discrete logarithms, quantum phase esti- 
mation, etc. 

Traditionally, QFT on r qubits is implemented by a 
quantum circuit consisting of 0(r) Hadamard gates and 
0(r 2 ) controlled-phase gates Q- Experimental demon- 
strations include synthesis of 3-qubit QFT in nuclear 
magnetic resonance (NMR) systems Q, order finding 
with NMR phase estimation with NMR [n|, Shor's 
factorization in NMR [ll|, in ion traps and using 
a "compiled version" of Shor's algorithm with photonic 
qubits [HI, [H[ • Further theoretical proposals for imple- 
mentations of QFT include atoms in cavity QED [15l ]. 
entangled multilevel atoms [l6| , trapped ions with House- 
holder reflections Hill, linear optics with Cooley- 
waveguide arrays [2(|, etc. 
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Tukey's algorithm 

The largest numbers factorized experimentally by 
Shor's algorithm hitherto are 15 [H OJ] and 21 [H[. The 
primary obstacle for demonstration of Shor's factoriza- 
tion for larger numbers is the large number of one- and 
two-qubit gates required. A "general-purpose" Shor's al- 
gorithm for an L-bit number demands L? gates and 5L+1 
qubits [22| ; an implementation using a linear ion trap 
would require about 396L 3 laser pulses [2^. "Special- 
purpose" algorithms that exploit special properties of the 
input number are much faster: the number 15 can be fac- 
tored with 6 qubits and 38 pulses only (22j. 

Another practical difficulty of the QFT algorithm is 
the use of two-qubit control-phase gates, which, for large 
number of qubits, involve very small phases. To this end, 
an "approximate" QFT has been proposed 0, [U, HH , in 
which the phase shift gates requiring highest precision 
are omitted. 

Griffiths and Niu proposed a "semiclassical" QFT, 
wherein the costly two-qubit gates are replaced by se- 
rial single-qubit rotations supplemented with classical 



measurements [23j. Such a semiclassical QFT has been 
demonstrated recently with three trapped ions [l2T |. 

In the present work, we propose to construct QFT by 
a novel approach which uses a special class of Hamil- 
tonians, having a circulant symmetry. Such Hamiltoni- 
ans have the advantage that their eigenvectors are the 
columns of the QFT (hence the latter diagonalizes the 
Hamiltonian), and they do not depend on the particu- 
lar elements of the Hamiltonian, as far as the circulant 
symmetry is conserved. This important feature allows 
one to construct QFT in a single interaction step; it also 
makes this techniques robust against variations in the in- 
teraction parameters. The present paper uses a similar 
approach as Unanyan et al. [26| . who proposed to use 
circulant Hamiltonians in order to create coherent super- 
positions of states. 



II. BACKGROUND 

Quantum Fourier transform. The A^-dimensional 
QFT is defined with its action on an orthonormal ba- 
sis |0),|l),...,|iV-l}: 
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It transforms a single state into an equal superposition 
of states with specific phase factors. The inverse QFT is 
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In a matrix form QFT is a square matrix with elements 
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Circulant matrix. An N x N matrix C of the form 



C = 
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is called a circulant matrix. It is a special case of a 
Toeplitz matrix [27j and it is completely denned by its 
first vector-column (or row). The other columns (rows) 
are just cyclic permutations of it. The circulant matrices 
have some very interesting properties. The most impor- 
tant one in the present context is that the eigenvectors of 
a circulant matrix of a given size are the vector-columns 
of the discrete Fourier transform © of the same size; 
hence they do not depend on the elements of the circu- 
lant matrix. The eigenvalues A„ of the circulant matrix, 
though, are phased sums of the matrix elements: 
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III. DESIGN OF THE HAMILTONIAN 

In order to synthesize QFT, we use a special time- 
dependent Hamiltonian of the form [26[ 



H(t) = /(t)Ho + ff(t)Hi, 



(6) 



where f(t) and g(t) are (generally pulse-shaped) real- 
valued functions, such that f(t) precedes g(t) in time, 
i.e. 







/(*) 



For instance, we can take 



f(t) = [1 - tanh(i/T)]/2, 
ff (t) = [l + tanh(i/T)]/2. 



(7) 



(8a) 
(8b) 



Therefore, the Hamiltonian © has the asymptotics 

H ~+^* H(t) *-=^ Hi. (9) 

We demand Ho to be a diagonal matrix in which the 
energies of all states (the diagonal elements) are non- 
degenerate 



U = diag(E 1 ,E 2 ,...,E N ). 



(10) 



For Hi we choose a circulant matrix, with the condition 
that the eigenvalues should be well separated from each 
other. Because the Hamiltonian has to be Hcrmitian, Hi 
is not a most general circulant matrix, but a Hermitian 
circulant matrix. 

Because the Hamiltonian © at t — > oo has a circu- 
lant symmetry, its eigenvectors are the vector-columns of 
QFT. However, each eigenvector |n) may have an adia- 
batic phase factor e IQ ™, acquired in the end of the evo- 
lution, which may be different for each \n). This means 
that for such a Hamiltonian ((6|) , adiabatic evolution will 
perform the QFT |T]) (possibly after renumbering of the 
basis states), but with some additional phases a n , 
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FIG. 1: Eigenvalues of the Hamiltonian (|6]) as a function of 
time for Ho and Hi given by Eqs. (|14a|) and ()14b|) . with 
V — E(l + i/3), whereas g(t) and f(t) are given by Eqs. 
dT6al) and flgbl. 



The phases a n are just integrals over the quasienergies, 
as follows from the adiabatic theorem [28|. The inverse 
Fourier transform would be 



JV-l 



(F^W) = --L J2 e~ iak e- 2vmk ' N \k), 



(12) 



fe=0 



and it can be accomplished by adiabatic evolution with 
the Hamiltonian 



H(t) = g(t)H + /(t)Hi. 
For example, for iV = 4we can have 

H = diag(-£, -E/3, E/3, E) 



Hi 



V V* 

V* V 

V* V 

V V* 



(13) 
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For laser-driven atomic and molecular transitions, the in- 
teraction V is given by the Rabi frequency fi: V = |M1 
Insofar as the eigenvalues of the circulant matrix are 
given by Eq. ([5]) , one has to choose the interaction energy 
V in such a way that the eigenenergies have well sepa- 
rated values. Another requirement for adiabatic evolu- 
tion is that the functions f(t) and g(t) change sufficiently 
slowly, so that the nonadiabatic coupling (xm(t) \ Xn(t)} 
between each pair of adiabatic states \Xm(t)} and \Xn(t)) 
remains negligibly small compared to the separation of 
the eigenenergies e m (t) and s n {t), 



,(t)-£„(t)| » \(Xm(t) | X n{t))\ 



1 

T' 



(15) 



where T is the interaction duration. 

In the numeric examples we use a hyperbolic-secant 
mask for the functions f(t) and g(t): 



fit) 

g(t) 



sech(f/r)[l 
sech(t/i-)[l 



-tanh(i/T)], 
■tanh(i/T)]. 



(16a) 
(16b) 
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These factors are chosen for implementation feasibility; 
they do not change the (all-important) asymptotic be- 
haviour of the eigenstates \Xn(t)). Figure [T] shows the 
evolution of the eigenvalues of the Hamiltonian ^ for 
V = i£(l + i/3). For this choice of V the eigenenergies are 
non-degenerate (except at infinite times, which is irrele- 
vant because there is no interaction) and the adiabatic 
evolution is enabled. 



IV. QUANTUM PHASE ESTIMATION 

We shall show now that the QFT propagator (fTTj) . 
which results from the Hamiltonian jB]), can be used to 
realize quantum algorithms, despite the presence of the 
adiabatic phase factors e la ™. We consider the quantum 
phase estimation algorithm [l[ , which is the key for many 
other algorithms, such as Shor's factorization. We briefly 
summarize here the essence of this algorithm. 

Let us consider a unitary operator U, which has 
an eigenvector \u) and a corresponding eigenvalue 
cxp(27ri</>), where <j> £ [0, 1). We assume that we are able 
to prepare state \u) and to perform the controlled-U 2J 
operation, for non-negative integer j. The goal of the 
algorithm is to estimate <p. To this end, the algorithm 
uses two registers. The first register contains r qubits 
initially in state |0) and the second one starts in state 
I it), containing as many qubits as needed to store \u). 

The procedure starts with the application of a 
Hadamard transform [l[ to the first register, followed by 
the application of controlled-U operations on the second 
register, with U raised to successive powers of two. The 
final state of the first register is 

2~ r / 2 (|0) + e 2m2r "^|l))i(|0) + c 2m2 ^ 2 ^|l)) 2 ■ • ■ 

2 r -l 

x(|0) + e 27ri2 °*|l)) r = 2-' r / 2 e 27nfc0 |fc), (17) 

fc=0 

and the second register stays in state \u). Now let us 
suppose that <f> can be expressed using a r-bit expansion 

= 0.^02...^ = ^ + ^ + --- + ^, (18) 

where O.</>i02 ■ ■ • 4>r represents a binary fraction. Then 
state (|17p can be written as 

2- r / 2 (|0) + e 27rl0 '^|l))i(|0) + c 27ri0 -^- 1 ^|l)) 2 ■ • ■ 

x(|0) + e a>no.fc*r"*r|i)) r . (19) 

Finally, wc apply the inverse QFT in order to obtain the 
product state \4>\ ■ ■ ■ <p r ) . In our case we apply the phased 
inverse QFT (TTJ) and find 

e- 10, W|&. ■■&.), (20) 

where a{<p) is an adiabatic phase that depends on <j). 
Since this global phase a has no physical meaning, a 
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FIG. 2: Upper frame: Field functions f(t) and g(t), Eqs. (|16a[) 
and (|16b|) . Lower frame: Fidelity of phase estimation during 
the inverse Fourier transform. The parameters are: <j> = 0.75, 
t — T, V = E(l + i/3), E = 10/T. 

measurement in the computational basis would give us 
exactly <f>. We note that if <f> cannot be written as a r bit 
expansion (|18p . this procedure can still produce a good 
approximation to <f> with high probability [l[ . 

In Fig. [2] we plot the probability of state ([20)) during 
the inverse Fourier transformation (|12p . This probabil- 
ity is evaluated by solving numerically the Schrodingcr 
equation for the Hamiltonian (|13[) and is used as a mea- 
sure of the fidelity. The figure shows that when the phase 
4> has an exact expansion as a binary fraction, the final 
probability tends to unity. 

V. IMPLEMENTATIONS 

In this section we discuss a few simple systems, which 
can be used for implementing the Hamiltonian (|6|). 

J' = i <-> J" = i system. As a first example we 
consider the system formed of the magnetic sublevels in 
a J' = i «-> J" = i transition shown in Fig. [3^a), 
where J is the total angular momentum of each level. 
We apply two linearly polarized fields (the two polariza- 
tion directions being perpendicular), the second one seen 
as a superposition of two circularly polarized fields (ct+ 
and cr_). By ordering the magnetic sublevels in the se- 
quence \m' = — h), \m" = |), \m' = h), \m" = —h), and 
by suitably tuning the strengths and the relative phase 
of the two independent fields, we can adjust the interac- 
tion elements of the Hamiltonian and produce the desired 
circulant form (|14b[) . We note that because of the dif- 
ferent signs of some of the Clebsch-Gordan coefficients, 
one should redefine one of the probability amplitudes by 
changing its sign. 

As we want to implement the full Hamiltonian ([6]), 
we also need to realize its first part Ho. It is especially 
important to remove the degeneracies between the mag- 
netic sublevels. This can be accomplished by using a 
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./' = 




J"' = 



m" = — 1 



J" = \ 




m" = 1 



J' = f * 



m' = 



/" = 1 



(c) " 




1 m ' = m ' 



FIG. 3: Systems, which can be used in order to realize a 
circulant Hamiltonian: (a) J' = | <-> J" = | system, (b) 
J' = J" = 1 <-> J'" = system, and (c) J' = 1 <-» J" = 1 
system, where the m' — <-» m" = transition is dipole 
forbidden. 



static magnetic field, which induces m-dcpcndcnt Zee- 
man shifts, and a far-off resonant laser pulse, which will 
cause Stark shifts. Let the energy splitting due to Zee- 
man shift be Ez (the same for both ground and excited 
levels). The Stark shifts are generally different for the 
two levels: E gt s and -E e ,s, where 'g' and 'e' stand, respec- 
tively, for ground and excited. Hence, in order to realize 
the Hamiltonian (|14ap . we need to solve the following 
algebraic system 



\E Z 
\E z 
-\Ez 

nE7, 



E g,S 
E g,S 

E e ,s 
E e ,s 



-E, 
-E/3, 
E/3, 
E, 



(21a) 
(21b) 
(21c) 
(21d) 



which gives Ez = E e< s = —E gt s = \E. Moreover, 
because the energies of Ho need not be exactly evenly 
separated, our method is robust against fluctuations in 
the field parameters. Making a reference to Eq. © and 
Fig. [2j we conclude that the Stark and Zeeman fields, 
with the time dependence f(t), have to be applied before 
the polarized laser fields, with time dependence g(t). 

J' = <-> J" = 1 <-> J'" = system. Another system 
with N = 4 is the diamond system depicted in Fig. [31(b). 
Here again two linearly polarized laser fields are needed, 
but now they have parallel polarization directions. One 
advantage of this system is that only magnetic fields are 
sufficient to realize the first part of the Hamiltonian. The 
disadvantage is that the two independent fields generally 
come from two different lasers, because of the different 
frequencies of the transitions. 

J' = 1 J" = 1 system. The J' = 1 <-> J" = 1 
system, depicted in Fig. [DJc), contains six coupled m 
sublcvcls. In this system the circulant symmetry occurs 



because m! = <-> m" = is a dipole forbidden transi- 
tion. By ordering the magnetic sublevels in the sequence 
\m' = -1), \m" = 0), \m' = 1), \m" = 1), \m' = 0), 
\m" = — 1) we obtain a Hamiltonian of the type 
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where Oi and O2 are the Rabi frequencies between, re- 
spectively, states with different m and states with the 
same m. Each Clcbsch-Gordan coefficient is incorporated 
in the respective Rabi frequency. The Rabi frequencies 
are complex (needed to avoid eigenvalue degeneracies), 
with a phase difference between the left and right cir- 
cularly polarized components. After a phase transforma- 
tion of the amplitudes, c n — > e 1/3 "c n , with suitably chosen 
phase factors we can make the Hamiltonian take the 
form of a circulant matrix. The selection of the phases 
)3 n amounts to solving a simple linear algebraic system. 

The first part of the Hamiltonian Ho can be realized 
with auxiliary magnetic and electric fields, as for the J' = 
\ <-> J" = B system. 



VI. CONCLUSIONS 

The intrinsic symmetry of circulant matrices allows 
one to design Hamiltonians that can produce a discrete 
Fourier transform on a set of quantum states in a natural 
manner and in a single step, without the need to apply 
a large number of consecutive quantum gates. The de- 
signed Hamiltonian has different asymptotics: it is a non- 
degenerate diagonal matrix in the beginning and a cir- 
culant matrix in the end (or vice versa); the time depen- 
dence that connects the two should be sufficiently slow in 
order to enable adiabatic evolution. The resulting uni- 
tary transformation, which this Hamiltonian produces, 
differs from the standard QFT by additional (adiabatic) 
phase factors in the matrix columns; we show, however, 
that one can still construct the quantum phase estimation 
algorithm, which is an essential subroutine in many quan- 
tum algorithms. We have presented examples of simple 
atomic systems, the Hamiltonians of which can be tai- 
lored to obtain circulant symmetry. The construction 
of large-scale systems with circulant symmetry requires 
the design of a closed-loop linkage pattern; for instance, a 
chain of nearest-neighbor interactions supplemented with 
a (direct or effective) interaction between the two ends 
of the chain. 
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